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Abstract 



We study the influence of finite kinematic boundaries on the induced gluon radiation from 

a fast quark in a finite size quark-gluon plasma. The calculations are carried out for fixed 

and running coupling constant. We find that for running coupling constant the kinematic 

^ , correction to the radiative energy loss is small for quark energy > 5 GeV. Our results 

^ [ differ both analytically and numerically from that obtained by the GLV group [6]. The 

O I effect of the kinematic cut-offs is considerably smaller than reported in [6] . 

O 

■^ ■ 1. It is very likely that parton energy loss due to the induced gluon radiation caused 

by multiple scattering [1, 2, 3, 4, 5, 6, 7] plays a major role in suppression of high-pT 

hadron spectra in heavy-ion collisions observed at RHIC [8, 9]. There is an attractive 

idea [10] to use this phenomenon (usually called "jet-quenching") to obtain information 

Qj- about the density of hot quark-gluon plasma (QGP) produced in AA-coUisions. Such 

^ ■ a jet tomographic analysis requires accurate methods for evaluating the induced gluon 

emission. In recent years this problem has been attacked from several directions. In [3] 

k> ! (see also [11, 12, 13]) we have developed a light-cone path integral (LCPI) approach to 

^ \ the induced radiation. The induced gluon spectrum was expressed through a solution 

of a two-dimensional Schrodinger equation in the impact parameter space with an imag- 
inary potential. This approach accounts for the Landau-Pomeranchuk-Migdal (LPM) 
effect [14, 15], finite-size and mass effects which are important for the QGP produced in 
y4y4-collisions. In [2, 5] the radiative energy loss was addressed using diagrammatic for- 
malism. Similarly to the LCPI approach the BDMPS approach [2, 5] expresses the gluon 
spectrum through a solution of a two-dimensional Schrodinger equation in the impact pa- 
rameter space. However, the BDMPS formalism applies only in the limit of strong LPM 
suppression. In this regime it is equivalent [5, 13] to the approach [3]. The GLV group 
[6] has developed within the soft gluon approximation an approach in momentum space 
which applies to thin plasmas when the mean number, N, of jet scatterings is small, and 
performed calculations accounting for the N = 1, 2, 3 rescatterings. 

For applications of the formalisms [3, 5, 6] to the tomographic analysis of experimental 
data on AA-collisions it is important to understand the limits of applicability of these 
approaches. In the analyses [3, 5, 6] the QGP is modeled by a system of Debye screened 
color centers [1], and parton scattering is treated in the small-angle approximation. The 
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LCPI [3] and BDMPS [2, 5] approaches, formulated in the impact parameter space, in 
addition, imply that the integration over the transverse momenta can be extended up 
to infinity ignoring finite kinematic boundaries. The small-angle approximation for fast 
partons moving through the medium should work well for parton energy E ^ hd, where 
Hd is the Debye screening mass which plays the role of a natural infrared cut-off and 
energy scale for parton scatterings in the QGP. For RHIC and LHC conditions, where 
fi£, ~ 0.5 GeV [16], it means -E > 3 ^ 5 GeV. It is however not clear whether the 
approximation of static color centers neglecting recoil effects is adequate for i? ~ 5 GeV. 
For such energies the kinematic constraints on the momentum transfer q < qmax ~ VEEt^, 
where Eth ~ 3T is the typical thermal energy of quarks and gluons in the QGP, may be 
important. The GLV group [6] has reported that the kinematic cut-offs suppress greatly 
the parton energy loss, AE. Even at -E ~ IOOO/id for the leading N = 1 contribution 
to AE for a homogeneous QGP with thickness L ^ 5 fm the authors have found the 
suppression ~ 0.5, and for E ~ lO/iu they give the suppression ~ 0.16. 

The approaches [3, 5] become inapplicable when the kinematic bounds become impor- 
tant. The GLV formalism [6], which does not treat accurately parton scattering near the 
kinematic limit, also cannot be used for quantitative calculations in this regime. Thus, if 
the kinematic effect were as strong as found in [6], the available approaches to the induced 
radiation would be inapplicable even at LHC energies. For this reason the kinematic effect 
merits further investigations. In particular, it is clearly desirable to study the effect of 
the running coupling constant. The decrease of the coupling constant near the kinematic 
bounds should act as a natural cut-off of large parton transverse momenta and diminish 
the role of the recoil effects. Another remaining open question is related to the different 
cut-offs for the initial and final partons. The authors of [6] have used for scattering of the 
radiated gluon the same cut-off in the momentum transfer as that for the initial parton. 
However, for soft gluons with x <^ 1 (hereafter x is the gluon fractional momentum) the 
Qmax is Considerably smaller than for the initial parton. In the present paper we address 
the role of the kinematic cut-offs accounting for the running coupling constant and dif- 
ferent cut-offs for the initial and final partons. The analysis is performed for the N = 1 
scattering which dominates the induced spectrum for RHIC and LHC conditions. We find 
that although the difference in the initial and final state cut-offs changes the analytical 
form of the induced spectrum, numerically the effect is insignificant. For fixed coupling 
constant the kinematic corrections become important for i? < 10 -i- 20 GeV, and for run- 
ning one the kinematic effect is small even at -E ~ 5 GeV. We find that the kinematic 
effect is considerably smaller than found in [6]. 

2. We consider a fast quark with energy E produced aX z = (we choose the z-axis along 
the momentum of the quark) traversing a medium of thickness L, which eventually splits 
into a gluon and final quark with the energies xE and (1 — x)E respectively. We assume 
that parton energies are much larger than the thermal quasiparticle masses in the QGP. 
The N = 1 induced spectrum can be represented in the form 

dP /■ , , ,da^^{x,z) 
—— = / dzniz) 







where n{z) is the number density of the medium (the summation over the triplet (quark- 



antiquark) and octet (gluon) color states is implied on the right-hand side of (1)), and 
151 

dx 



— — ^^^ is the in-medium (z-dependent) Bethe-Heitler cross section. It can be written as 



] — Jbb + Jcc + Jdd + 2 Jfjc + 2Jcrf + 2Jdb 

dx 

+'^Jae + 2Ja/ + 2Jag + 2Jah (2) 

with Jq,^ given by 

E f 

Jap = -j^^^f^^y dqdpTo,{q, p)T^(q, p) , (3) 

where the amplitudes T^ diagrammatically are shown in Fig. 1, q and p are the transverse 
momenta of the t-channel and emitted gluons, respectively. Note that the interference 
between the double-gluon exchange diagrams (e), (f), (g), (h) and the diagram without 
gluon exchange (a) is important to ensure unitarity. 

The diagrams of Fig. 1 can be evaluated with the help of the ordinary perturbative 
formula _ 

T = j^ dz' J dpgi^fip, z')YA^{p, z')^,{p, z') , (4) 

where p is the transverse coordinate, ilJij{p,z') are the wave functions of the initial and 
final quarks, and A^{p,z') is the wave function of the emitted gluon (hereafter we omit 
the color factors and spin indices). In (4) we do not explicitly indicate the dependence 
of the wave functions on the position of the scattering center. The quark wave functions 
using the light-cone spinor basis can be written as 

^i(p> z') = exp{iEjz')Uj(j)j{p, z!) , (5) 

where the operator Uj reads 

Here Xj is the quark spinor (normalized to unity), a = 7^7, {3 = 7'^, and p = — iVj.. 
The transverse wave function (f)j{p, z') entering (5) is governed by the two-dimensional 
Schrodinger equation in which z' plays the role of time 
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d<P,{p,z') _ \{v^ + ml) 



dz' 
where 



2E, 



+ v{p,z 



(P, ^') , (7) 
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v{p, z') = 5{z' -z) j ——^ exp (iqp)t;(q) , t;(q) = - 



'^^^l 



is the potential generated by the one-gluon exchange between quark and the Debye 
screened color center. In the longitudinal direction we treat the potential as a point- 
like. In the same form one can represent the gluon wave function (up to an obvious 
change of the spin operator and color factors). 



The amplitudes entering (3) can be easily obtained from (4) treating in (7) the po- 
tential f as a perturbation. For the diagrams with gluon exchanges in the z' regions 
^ < z' < z and z' > z the transverse wave functions age given by the plane waves (with 
different transverse momenta in these two regions of z') 



[p, z ) oc exp 



PiP 



iPj + ^3 



IE, 



(9) 



The color center acts as a kick which changes the quark (or gluon) transverse momentum 
at z' = z. The corresponding amplitude oc f (q) for one-gluon exchange diagrams, and 
oc J dpv{p)v{q — p) for the double-gluon exchange ones. Note that, eventually, the p- 
integration in (4) ensures conservation of the transverse momentum. 

To account for the kinematic boundaries we introduce in the amplitudes the cut-off 
factors. For each t-channel gluon we modify the propagator introducing the cut-off factor 
d{<lm.ax — q)- Here Qmax IS the upper kinematic bound on the momentum transfer for 
the parton to which the t-channel gluon is attached. Also, we modify the qqg-vertex 
for splitting the initial fast quark into quark-gluon system introducing the cut-off factor 
0{kmax — k), k is the transverse momentum of the gluon in the frame where the total 
transverse momentum of the quark-gluon state equals zero, and kmax = E min(x, 1 — x) 
(here E is the initial quark energy). The above prescription ensures that parton scattering 
angles are small, and the momentum transfer does not exceed the kinematic bounds. 

Using (2)-(8) after straightforward but a bit cumbersome calculations the effective 
Bethe-Heitler cross section can be represented in the form 



da^^{x,z) daf"{x,z) 



dx 



dx 



+ 



da^^{x,z) 
dx 



(10) 



where 



daf^{x,z) 
dx 



F{k,q,z) 



tt'^x 



dddkj-^^^ [e{q, - q)F{k, q, z) 



+0{qi - q)F{k, q(l -x),z)- —Oiq^ - q)F{k, qx, z)] , 



m 



k'Q\k) (k-q)ke(k)e(k-q) 



(k' + e2)2 (k^ + e2)((k-q)2 + e2)_ 



1 — cos 



IZ 



l{k,x) 



(11) 
(12) 



da^^{x,z) 
dx 



2Ct f, ^x^ 
Ca7c^x[ "" 2 



dq , f^^'^l [CF{e{qo - g) - e{q, - q)) 



+CAieiq2-q)-9iqs-q))]- j dk 
with the following shorthands: 



q2 + ^ilf 



(k^ + e2)2 



/(k,x) 



2Ex{l -x) 
k2 + e2 



1 — cos 



IZ 



Jik,x) 



(13) 



(14) 



go = qmax{E), gi = qmax{Ex), q-2 = qmax{E{l - x)), gs = min(gi,g2), (15) 



0(k) = J as{k)9{kmax — k), e^ = m^x^ + m^(l — x), nig^g are the thermal quark and 
gluon quasiparticle masses, Ct,f,a are the color Casimir factors of the color center, quark 
and gluon respectively. Eq. (15) corresponds to the above described scheme when each 
scattered parton has its own g-cut-off factor. Note that in the soft gluon limit x <^ 1 
our formulas do not reduce to that of Ref. [6]. If one uses for the final partons the same 
Qmax as for the initial quark as was done in [6], the second term on the right-hand side 
of (10) vanishes. This term emerges inevitably because the initial and final partons have 
different phase space for their scattering. Below for comparison with [6] we also present 
the results for g, = QmaxiE) as in [6]. 

The quantity Lf = l{\<. = 0,x) characterizes the longitudinal scale of gluon emission, 
i.e., the gluon formation length. The induced spectrum depends crucially on the ratio 
Lf/L [3, 4, 17]. For gluons with small formation length Lf <^ L the finite-size effects 
caused by the oscillating cosine on the right-hand side of (12), (13) becomes small. In 
this regime the rapidly oscillating cosine can be neglected, and the effective cross section 
(10) becomes equal to the ordinary Bethe-Heitler one, i.e., to the cross section for a quark 
which approaches the color center from outside. In contrast, when Lj > L the finite size 
effects due to the cosine in (12), (13) suppress greatly the radiation rate as compared to the 
Bethe-Heitler one [4, 17]. This suppression, physically, is connected with small transverse 
size of the qg system (it is oc L). In this regime the t-channel gluons cannot distinguish 
the |g) and \qg) Fock components of the physical quark and for this reason the gluon 
emission turns out to be suppressed. One remark regarding the Bethe-Heitler regime for 
Lf <^ L is in order here. Diagrammatically, the ordinary Bethe-Heitler cross section is 
given by the diagrams (b), (c), (d) of Fig. 1 involving only one-gluon exchange. However, 
our formulas include the interference between the diagram (a) and (e), (f), (g), (h). The 
explanation of this fact is as follows. For a quark incident on the color center from outside 
the amplitudes (b), (c), (d) should be evaluated integrating over z' in (4) from — oo (with 
usual adiabatic switching off of the coupling constant for l-z'l ^ oo). For a quark produced 
in a hard reaction aX z' = the ^'-integration region is (0, oo). This gives rise to additional 
endpoint terms (corresponding to z' = 0) in the cross section which are absent when the 
lower limit equals — oo. However, similar endpoint terms emerge for the interference term 
involving the double-gluon exchange diagrams as well. They cancel exactly the endpoint 
terms stemming from the graph (b), (c), (d). As a result for z ^ oo our effective cross 
section (10) equals the ordinary Bethe-Heitler one. 

It should be noted that our method (and any other one based on the GW model [1] 
for the QGP) can only give an estimate for the kinematic correction. It is inapplicable in 
the regime of strong kinematic suppression when the spectrum becomes very sensitive to 
the detailed form of the kinematic cut-offs. This fact is closely connected with the anti- 
leading log character of the q, k-integrations in (11). Contrary to the ordinary leading 
log situation, say, in 7* -^ qq transition in deep inelastic scattering, where the typical 
values of the momentum transfer g is smaller than the internal momentum k, in the case 
of the induced gluon emission in the high energy limit when Lf ^ L the dominating 
contribution to the induced spectrum comes from g > /c. A detailed discussion of this 
phenomenon is given in [17]. 



The effective Bethe-Heitler cross section evaluated without kinematic cut-offs, i.e., 
with Qi = kmax = oo, is givcu by the first term on the right-hand side of (10). After 
the Fourier transform it can be represented in the impact parameter space in the form 
obtained previously [18, 17] within the LCPI formalism 

da^^ix z) f 

-^-^^ = Re j dp^*{p,x)a:i{p,x)'^m{p.x,z) , (16) 

where \E'(p, x) is the ordinary light-cone wave function for the q -^ gq transition in vacuum, 
\E'm,(p, a;, 2;) is the z-dependent light-cone wave function describing the quark-gluon Fock 
component at the longitudinal coordinate z, and a^{p,x) is the three-body cross section 
of a qqg system with a particle in the medium (the explicit form of the wave functions 
and three-body cross section can be found in [17]). In the qqg system antiquark is located 
(in the transverse space) in the center of mass of the qg pair, and the relative separations 
satisfies the relation (p^ — p^x = (1 — x){pq — p^). 

3. We have performed numerical calculations for fixed and running coupling constant. In 

the first case we take Ug = 0.5 [19]. For running coupling constant we use the one-loop 

formula with Aqcd = 0.3 GeV frozen at the value a^ = 0.65. In this case a^ approximately 

satisfies the relation 

r-2 GeV a (k) 

dk^^^-^ ~ 0.36 GeV (17) 

TT 

obtained from the analysis of the heavy quark energy loss [20]. We have carried out the 
calculations for expanding plasma. We use the Bjorken [21] model with Tt"^ = TqTq, 
and take the initial conditions suggested in [22] for heavy ion central collisions at RHIC: 
To = 446 MeV and tq = 0.147 fm. For the upper limit of the 2;-integration in (1) we take 
L = Ra ~ 6 fm. For quark and gluon quasiparticle masses we use the values obtained 
in [16] from the lattice data rriq ^ 0.3 and rrig ^ 0.4 GeV. With the above value of rUg 
from the perturbative relation jio = \plvag one obtains for the Debye screening mass 
/iz) ~ 0.57 GeV. For the energy dependence of the maximum momentum transfer we take 
(^^^^{E) PS EEth with Eth = 750 MeV. It is smaller than q^^^{E) ^ 3E/iD used in [6]. 

In Figs. 2, 3 we plot the induced gluon spectrum for i? =5, 10 and 20 GeV evaluated 
using (1), (10)-(15) with (solid line) and without (dashed line) the kinematic cut-offs 
for fixed and running q;^. The results with kinematic cut-offs have been obtained for qi 
given in (15). For this version we also plot the spectrum without the second term in (10) 
(dotted line). In Figs. 2, 3 we also show the results for the same kinematic cut-offs for 
initial and final partons obtained with qi = qmax{E) (long-dashed line). From Figs. 2, 3 
one sees that the kinematic cut-offs become especially important when the energy of the 
radiated gluon (or the final quark) < 1 ^ 2 GeV. The kinematic correction is smaller for 
running coupling constant. It is also seen that for fixed coupling constant the relative 
contribution from the second term in (10) is larger. It is natural since for the running ag 
the contribution from large transverse momenta is suppressed. The total spectrum in the 
above two scheme of the g-cut-off turns out to be approximately the same. 

To illustrate the effect of the kinematic cut-offs on the quark energy loss in Fig. 4 we 
plot the energy dependence of the kinematic K-iactor 

K{E) = ^, (18) 



where ^Ef^b, and AEi,b. are the quark energy losses evaluated with (for qi as given in 
(15)) and without kinematic constraints, respectively. We define AE as 

AE = E dxx— (19) 

with Xmin = ^TT'g/E, Xmax = l—mq/E. Fig. 4 demonstrates that for fixed coupling constant 
the kinematic cut-offs are important for £" < 20 GeV. For running coupling the kinematic 
correction is small even for i? ~ 5 GeV. One can see that the difference between the 
cut-offs given in (15) (thick curves) and q^ = q^axiE) (thin curves) is less than ~ 10%. 

The kinematic K-iactor shown in Fig. 4 is obtained for an expanding plasma with 
n{z) oc 1/z. For comparison with the analysis [6] we have also carried out calculations 
for a homogeneous plasma. In this case the kinematic effect is weaker since the relative 
contribution of the region of small z where the typical parton transverse momenta are large 
(they are oc l/^/i) is smaller. For correspondence with [6] we have taken g^^^ = SEfio 
(the same for the initial and final partons) and k'^^x — 4i5^ min(x^,a;(l — x)) used in 
[6] which give somewhat smaller kinematic effect than our cut-offs. We have obtained a 
small kinematic effect, say, K ^ 0.9 at E = 5 GeV and K ^ 0.94 at -E = 10 GeV. (we 
have used fixed coupling constant as in [6]). It is considerably larger than the suppression 
reported in [6] (~ 1/6 for E = 5 GeV). In connection with strong kinematic suppression 
reported in [6] we would like to emphasize one more time that within the approximation 
of static color centers [1] when the kinematic cut-offs are introduced by hand the regime 
of strong kinematic effect cannot be described accurately. It is clear that the analysis of 
the induced radiation in this regime requires an accurate treatment of the recoil effects. In 
this case the fast partons moving through QGP and partons from QGP should be treated 
on an even footing. Note also that in this regime suppression of the radiative energy loss 
may largely be compensated by the collisional energy loss due to strong recoil effects. 

In summary, the form of the induced gluon spectrum obtained in the present analysis 
shows that the kinematic effect is relatively small and is mainly important near to the 
endpoints a; ~ and a; ~ 1 when the energy of the radiated gluon (or of the final quark) is 
about ~ 2^3 units of the Debye screening mass, i.e., about 1-^2 GeV for RHIC and LHC 
conditions. For fixed coupling constant the kinematic correction to the quark energy loss 
becomes small for i? > 20 GeV, for running coupling constant it is small even at i? ~ 5 
GeV. The kinematic effect found in our analyses is considerably smaller than reported in 
[6]. Our results say that in the region of the gluon fractional momentum S <x <1 — 6 
{6 ~ (2-=-3)/i£)/i5), the induced spectrum can be evaluated to reasonable accuracy within 
the LCPI approach [3] which ignores the kinematic bounds. This approach can be used 
for evaluation of the energy loss and nuclear suppression factor [19] for RHIC and LHC 
energies. 
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Figure 1: The set of the Feynman diagrams for the A^ = 1 contribution to the induced 
gluon spectrum. 
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Figure 2: The spectrum of the induced q ^ gq transition versus the gluon energy xE 
for RHIC conditions for fixed couphng constant. The sohd hues are for qi given in (15) 
and the long-dashed hues are for the same g-cut-offs for the initial and final partons with 
qi = qmax{E). The dashed lines show the spectrum obtained without kinematics cut-offs. 
The spectrum without the second term in (10) for qi given in (15) is shown by the dotted 
curves. 
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Figure 3: The same as in Fig. 2 but for running coupling constant. 
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Figure 4: The energy dependence of the kinematic K-iactor (18) for RHIC conditions for 
running (sohd hne) and fixed (dashed hne) couphng constant. The thick hues are for the 
g-cut-off given in (15) and the thin hues are for qi = qmax{E). 
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